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Abstract. We make a detailed study of the first and second-order SUSY partners of a one- 
dimensional free Hamiltonian with a singular perturbation proportional to a Dirac delta 
function. It is shown that the second-order transformations increase the spectral manipula- 
tion possibilities offered by the standard first-order supersymmetric quantum mechanics. 
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1 Introduction 

The study of one-dimensional Hamiltonians with a point interaction has received renewed atten- 
tion during the past two decades [Il[2l[3l[ll[5l[6l[7j- In general, a point interaction is described 
by a potential concentrated either in a single or a discrete number of points as it happens, e.g., 
for the Dirac delta or its derivative. Mathematically, in order to define these potentials, we use 
the theory of extensions of symmetric operators with equal deficiency indices. These extensions 
have domains which are characterized by some matching conditions for the wave functions at 
the points supporting the interaction [21 [HI O [lOl [11] . In particular, the Dirac delta barrier or 
well have been extensively studied in this way with or without other interactions ^12j, with or 
without mass discontinuities at the singular points etc. [131 [Hj. 

On the other hand, supersymmetric quantum mechanics (SUSY QM) has emerged as the 
standard technique for generating new potentials with known spectra departing from an initial 

one[I3[M[IZl[M[Il[2ni[2Il[23[231[H[23[2a[2Zl[251[2a[3a 

[39l HOl E] . The method has been applied successfully to regular one-dimensional potentials 
defined on the full real line [421 EH] , on the positive semi-axis [44^ I45| or in a finite interval [l6] . 
Although there are some works dealing with SUSY QM applied to point potentials [171 UHl [iHl 
[50I [STl [52I [53] . however the corresponding study has been done just for particular first-order 
SUSY transformations, without analyzing the full possibilities of spectral manipulation offered 
by the method. It is interesting to note as well that a point potential may appear as hidden 
super symmetries [5^ [55]. 

Now, it is the appropriate time for studying the behavior of point potentials with bound 
states under SUSY QM. Due to the calculation complexity, we shall focus our attention to first 
and second-order transformations, which anyway are interesting by themselves [56 [ I57 [ \58 \ 159] . 



*This paper is a contribution to the Proceedings of the Workshop "Supersymmetric Quantum Me- 
chanics and Spectral Design" (July 18-30, 2010, Bcnasque, Spain). The full collection is available at 
|http://www.emis.de/journals/SlGMA/SUSYQM2010.html| 
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We shall restrict the discussion to the following one-dimensional Hamiltonian 
1 £ 

^0 = -^-^ + yo{x), Vo{x) := -a6{x), a > 0, (1) 

which is mathematically well defined and self-adjoint provided that we use as its domain T> the 
subspace of the Sobolev space PV2^(]R/{0}) such that for any ip{x) G V, one has: 

(2) 

where 'ijj{0+), tp'{0+) and -0(0—), 'ip'{0—) are the right and left limits of ip{x), 'ip'{x) at the origin 
respectively [2]. 

In order to achieve our goal, we have organized this paper as follows: in Section [21 we will 
study the solutions of the stationary Schrodinger equation for the Hamiltonian Hq given by ([1]). 
In Section [3] we will apply the first-order SUSY techniques, in Section H] we will analyse the 
second-order transformations and in Section [5] we will present our conclusions. 

2 Solution of the Schrodinger equation 

Let us evaluate in the first place the general solution of the stationary Schrodinger equation for 
an arbitrary e = — A:^/2 < 0: 

Hqu{x) = eu{x), (3) 

with Hq given in ([T]). There is one solution vanishing for x — ?• — oo, denoted of the form 

u+{x) = e^''H{-x) + (a e^"" + f3e-''^)H{x), k > 0, (4) 

where H{x) is the Heaviside step function, and a, (3 are constants to be determined from the 
discontinuity equations ([2]) . We need as well the derivative of (x) , 

nV(x) = ku+{x) - 2kPe-^''H{x) + (a + /3 - l)6{x). (5) 
From equations ([4]) and ([5]) it turns out that 

n+(0-K) = a + /3, n+(0-) = 1, u'j^{Q+) = k{a - (3) , v!_^{Q-) = k. (6) 
On the other hand, using equations ([2]) and ([6]), we obtain: 

Q-h/3 = l, a- 13 = 1 - 2a, 
where a = a/k. Hence: 

a = 1 — a, j3 = a. 

Inserting these expressions in equations ^ and ([5]), we finally get 

u+{x) = e^'' -h{e^'' -e-^'')H{x), (7) 
u^{x) = ku+{x) - 2ae~''^H{x). 

Note that the Hamiltonian Hq in equation ([T|) is invariant under the change x — )• —x. Thus, 
we can find a second linearly independent solution n„(x) for the same e = —k'^/2, vanishing 
now for X — 7- oo, by applying this transformation to u^{x): 

= a(e'=^-e-'=^)i/(-x) + e-'=^. (8) 
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Moreover: 

u'_ (x) 



-A:n_(x) + 2ae*^^F(-x) 



Finally, the general solution of equation ([3]) for e = —k'^/2 < is a linear combination of both d?]) 
and dHI) which, up to an unessential constant factor, becomes: 



(9) 



u{x) = u+{x) + Du.{x) = e'^^ + L'e-'^^ - S(e^^ - e-^^^j [H{x) - DH{-x)] , 
where D is a constant. The corresponding derivative is given by: 

u'{x) = -ku{x) + 2A;e''^ [1 - dH{x) + DaH{-x)] . (10) 
Note that, up to normalization, both solutions u±{x) lead to the same bound state for ko = a: 

Mx) = V^[e'""'H{-x) + e-'"'''H{x)] . 
The corresponding eigenvalue becomes 



En 



a 



which coincides with the result derived in 

On the other hand, the scattering states for e = /2 > can be simply obtained from 
the solutions given in equations (O, (HI) by the substitution k — ?• —in, k > 0. In particular, 
for a probability flux approaching the singularity from — oo the corresponding scattering state 
arises in this way from the u_(x) of equation which (up to unessential constant factor) leads 
to 



ip{x) 



la 



K — la 



H{-x) + 



K — la 



-e"'''H{x). 



(11) 



It is clear now that the reflection R and transmition T coefficients become the standard ones 
(see, e.g., [60]): 



R 



K — la 



T 



K — la 



+ d 



2 • 



(12) 



3 First-order SUSY transformations 



Let us start with the initial Schrodinger Hamiltonian Hq given in ([T]). As it is well known (see, 
e.g., [2Z1 HI] and the references cited there), its first-order SUSY partner, 

is intertwined with Hq in the way 

HiA+ = A+Ho, (13) 

where 
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Here, the transformation function u{x) is the seed sokition given in associated to the fac- 
torization energy e = —k'^/2 and satisfying equation ([3]). The SUSY partner potential Vi{x) of 
Vo(x) is given by: 



Vi{x) = Vo{x) - [lnu{x)]" . 



(14) 



We assume the standard restriction e < i^o =^ ^' ^ ^o, in order to avoid the creation of new 
singularities in Vi{x) with respect to those of Vq{x). Note that, from equation (|14p and ([3]) we 
have: 



Viix) = Voix) 



u"{x) 
u{x) 



+ 



u'{x) 
u{x) 



-Vo(x) + 2e + 



Hence, a straightforward calculation using equations 



u'{x) 
u{x) 

lUD leads to: 



u'{x) 
u{x) 



4A:2(1 



v?{x) 



■[D + aD'^H{-x) + aH{x)\. 



As e = —k'^/2, equations ([T5]), ([T6|) give: 



Vi{x) =a5{x) 



4A:2(1 - a)[D + aD^H{-x) + aH{x)] 
{e^^ + De-'''' + 2a smh{kx)[DH{-x) - H{x)]y ' 



(15) 



(16) 



(17) 



Note that the denominator of equation (jl7p never vanishes for x E (—00,00) and D > 0. 
Moreover, it can be seen that the delta term in Vi{x) is now repulsive (since a > 0). 

A straightforward consequence of the intertwining relationship (jl3p is that for any eigenfunc- 
tion of Hq associated to the eigenvalue E {HqiI) = Eip) such that A^if) / 0, it turns out that 
oc A^iIj oc W {u,^) /u is a corresponding eigenfunction of Hi associated to E. Moreover, 
if il: satisfies as well equation ([2]) it turns out that ij)^^^ now obeys: 

( ^^(^^0+) \ 
1^ V^^')'(o+) I " 

which is consistent with the fact that the intensity of the delta term in Vi{x) has an opposite 
sign compared with Vq{x) and the second term of Vi{x) has just a finite discontinuity at x = 
(see equation p!7|) ). 

Concerning the spectrum of Hi , let us note in the first place that transforms the scattering 
eigenfunctions of Hq into the corresponding ones of Hi. In particular, the wavefunction tjj{x) 
given in equation (jllh . when transformed by acting on it with , produces an expression 
ip^^\x) which is a bit large to be presented here. However, for large values of that expression 
reduces to the following scattering one (up to a constant factor): 




il)^^\x) 



\x\- 



a 



+ 



-k + iK 
a + IK J \ k + IK 
IK \ I —k + IK 



H{-x) 



a + IK 



k + IK 



H{x). 



This means that the initial reflection and transmission coefficients are unchanged under the 
first-order SUSY transformation (compare equation (jl2p ). We thus conclude that the continuous 
spectrum of Hq belongs as well to the spectrum of Hi. 

Let us note that the differences in the spectra of Hi and Hq rely in general in the modifications 
produced by a non-singular SUSY transformation on the discrete part of the initial spectrum. For 
first-order transformations, these changes can be classified according to the essentially different 
combinations of the parameters > and k > kQ which characterize the seed eigenfunction u{x). 
We can find three different situations. 
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(i) Creation of a new ground state at e < Eq. This case appears for D > 0, k > kQ. Here, the 
eigenfunction ■ipi^'^ oc l/u{x) of Hi associated to e is square-integrable. Moreover, since the 
mapped initial ground state tp^^^ = ^ W{u^a) a normalized eigenfunction 
of Hx with eigenvalue Eq, then Sp(i/i) = {e, E^} U [0, oo) = {e} U Sp(Fo)- 

(ii) Isospectral transformations. These are achieved from the previous case either by taking 
D — 7- or D — 7- oo. Since in both situations u{x) goes to zero at one of the ends of the 
x-domain, it turns out that ipi^^ oc l/u{x) is no longer square-integrable, although tp^^ is. 
Thus, Sp(Fi) = {^o} U [0,oo) = Sp(Fo). 

(iii) Deleting Eq. This situation arises from the previous one by taking k = kQ = a [a = a/k = l). 
Since u{x) oc iIjq{x) is square-integrable, then ip^^ oc \/u{x) is not normalizable, and then 
Sp(//i) = [0, oo). From equation (fT7|l . it is clear that now 

Vi{x) = a6{x). (18) 

This means that, by deleting the bound state of the attractive delta well Vb(x) = —a5{x), 
a > 0, which is placed at Eq = — we recover the repulsive delta barrier of equa- 
tion (jlSp . a standard result well known in the literature. 



4 Second-order SUSY transformations 

In this section it will be illustrated, by means of the delta-well potential, the advantages for 
manipulating spectra of the second-order SUSY transformations [561 1571 [58l |59] compared with 
the first-order ones. It is nowadays known that the second-order SUSY partners H2 of the initial 
Hamiltonian Hq can be generated either by employing two eigenfunctions tti(x), U2{x) of Hq, 
not necessarily physical, associated to two different factorization energies £1^25 ci 7^ ^2 [371 |4T] 
or by an appropriate eigenfunction ui{x) in the limit when 62 — ?• ei (the so called confluent 
case |61l I62j). In both situations the two Hamiltonians Hq, H2 are intertwined by a second- 
order operator in the way 

H2B^ = B^Hq, 

where 




ui{x) ' 



^2 

the new Hamiltonian H2 takes the standard Schrodinger form 

and the second-order SUSY partner V2{x) of the initial potential Vo(x) is given by 

V2{x) = VQ{x)-[lnw{x)]", (19) 

the real function w{x) being proportional in general to the Wronskian of two generalized eigen- 
functions of Hq |63j . An explicit classification of the several second-order SUSY transformations 
is next given. 
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4.1 Confluent case HH [62] 

Let us consider in the first place the limit £2 — ei = e = —k'^/2 < 0, taking as seed the 

Schrodinger solution u+{x) vanishing as x 1 — > —00, which means to take the u{x) given in 
equation Q with D = 0, namely: 

u{x) = e^^ - 2a smh{kx) H (x) . (20) 

In this case the real function w^x) appearing in equation (fT9|) takes the form [62] 

v?{y) dy. 

'OO 

An explicit calculation for x < leads to: 

[■X ^2kx 

00 



u\y)dy= / e'''ydy = — 



On the other hand, for x > it turns out that: 

r /"O rx 

u'^{v)dy=l u^{y)dy+ j u^{y)dy 
00 J —00 Jo 

a (1 - 2d) 2fcx , «^ 



+ ^ — ^ e^''^ + — sinh(2A;x) + 2a(l - d)x. 



By combining these two results, we obtain 



w{x) = Wo + + 



a a 2kx a 



■ H{x). (21) 



^2kx _^ sinh(2fcx) + 2a(l - d)x 
k k k 



The second-order SUSY partner potential of Vb(x) becomes now 

. X N , ^"^(3;) 2u{x)u'{x) 
V2{x) = -a5{x) + . (22 

Note that, since n'(0+) = k — 2a and n'(0— ) = /c, then ti'(x) and consequently the potential 
difference AV{x) = V2(x) — Vb(x) have a finite discontinuity at x = 0. 

In order to avoid the arising of extra singularities for V2(x) with respect to Vo(x) we have 
to take vjQ > 0. Concerning the spectrum of H2, a similar calculation as in the first-order case 
shows that the scattering eigenfunctions of Hq are mapped into the corresponding ones of H2, 
i.e., the energy interval [0, 00) belongs to Sp{H2). As for the discrete part of the spectrum, 
several possibilities of spectral manipulation emerge, according to how we choose k and wq. 

(i) Creating a new bound state at e ^ Eq. This case appears by taking wq > and k ^ k^ = a. 
Since 

lim ■(/'e^^ lim — 7~T e"'"'^', 

|x|-s>oo |x|-s>oo w{x) 

(2) 

the eigenfunction ip^ of H2 associated to e is square- integrable, i.e., a new bound state 
has been created at e, either below the ground state for fc > fco or above it for k < k^. 
The last option is illustrated in Fig. [H where we have plotted the potential difference 
AV{x) as a function of x for a = 2, k = 1, wq = 1, i.e., a new level was created at 
e = — l/2>£'o = — 2 (see the two gray horizontal lines in the same graph). Note the 
existence of a finite discontinuity in AV{x) at x = 0, induced by a similar discontinuity 
of u'{x) at the same point. 
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3 




Figure 1. Potential difference AV{x) as function of x (blue lines) induced by the confluent second order 
SUSY transformation for a = 2, fc = 1, wq = 1- A new level was created at e = —1/2, above the initial 
ground state Eq — ~2 (gray horizontal lines). 

(ii) Isospectral transformations. They arise in the first place as a limit of the previous case 
for e ^ Eq and wq — t- 0. Note that the long explicit expression for the V2{x) of (I22p which 
would appear if we would substitute explicitly the u{x) and w{x) of equations (j20p . (j2ip 
becomes strongly simplified in this limit: 

V2{x) = -a5{x) + H{x) {Sdk^e^^''[{a - l)e^^'' - d]{{d - l)[2kx{a - 1) + 1 - 2d]e^^^ 
+ d[2kx{d - 1) + 25 - 3]}) /{(a - l)2e^^^+ 25[1 - 2kx{d - l)]e^'''' - d^f . 

Since now 

u{x) 



Hm 

x^-co w[x) 



OO, 



it turns out that e Sp(i?2) = {£^0} U [0,oo) = Sp{Ho). 

An alternative way to produce isospectral transformations is to use the single bound state 
of i^o for evaluating w{x). The corresponding formula is achieved from equation (|2ip by 
taking k = kQ = a, d = 1, which leads to: 

„2ax 



w{x) = wq -\ sinh^{ax)H(x) . 

2a a 



Hence 



V2{x) = —a5{x) 



TylH{-x) + 



{2woa + e2«^')^ ^ ' [2(1 + t(;oo)e2«^ - If 
Note that now w{x) does not have any node for 
1 



H{x). 



(23) 



(24) 



wq G 1^—00, J U (0, 00) . 

Moreover, in this domain it turns out that 
hm 44oce-N, 

i.e., V'e u{x)/w{x) is square-integrable Sp(i/2) = {^o} U [0, 00) = Sp(i7o) 
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(iii) Deleting the ground state of Hq . By taking now the limit of equation (j23p for itjQ — )• or 
wq —7- —1/a, it turns out that Um u{x)/w{x) = oo or lim u{x)/w{x) = oo respectively. 

x^—oo I— >-oo 

(2) 

In both cases ■0e is not square-integrable and then 
^0 0Sp(i72) = [0,oo). 

This result means that we have deleted the ground state of Hq in order to obtain H2. For 
1/70 —S" the potential of equation (p^ becomes 

V,{x) = -a5{x) + ^^^2ax_^y Hix). (25) 

On the other hand, for wq — )• —1/a the corresponding potential V2(x) is obtained from the 
previous one by the change x — )• —x. 



Let us remark that, although the final spectra of the SUSY partner Hamiltonians of Hq 
are the same when deleting its ground state in the first-order and in the confluent second- 
order transformations, however the potentials Vi{x) and V2{x) are physically different (compare 
equations (fTSj) and (j25]) ). In particular, note the opposite signs of the coefficients of the Dirac 
delta function for both potentials. 



4.2 Complex case |64ll65ll66] 

Let us assume that k = kn + ikj is complex with kn > 0, kj G M, and suppose that the two 
involved factorization energies are now given by e = —k'^/2 and e, where z denotes the complex 
conjugate of z. Since we need to avoid the arising of extra singularities in the new potential, we 
will take a Schrodinger seed solution vanishing at one of the ends of the x-domain in the form 
given in equation (I20p with k £ C, namely, 

u{x) = e''^ - ^ smh{kx)H{x), u{x) = - sinh{kx)H{x). (26) 
k k 

To compute now the second-order SUSY partner potential V2{x), we have to obtain in the first 
place the Wronskian W{u, u) and then the real function 

W{u,u) 

This calculation is cumbersome but otherwise straightforward, which leads to: 

^^^^ ^ ~2k^ ^ \ ~ 1^ [^°^^(^^«^) ~ cos(2A;/x)] 

+ ]%~k^kj sinh(2fcfix) - kRsm{2kix)] ^H{x). (27) 

Then, V2{x) will be given by 

\u{x)\'^ [u{x)u' (x) + u{x)u' (x)] 
V2{x) = -ad{x) + , (28) 

iW^^Xj w[x) 



with u{x), u{x) and w{x) as given in equations (|26p and (|27p . An illustration of the potential 
difference AV{x) as function of x for a = 2, A; = 1/100 + i/10 is given in Fig. [2j 
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Figure 2. Potential difference AV^(a;) as function of x induced by the complex second order SUSY 
transformation for a = 2, fc = 1/100 + i/10. The two potentials V2{x) and Vo(x) are isospcctral. 



Note that these equations become highly simplified if = a: 

w{x) 



2a V CL 



° ,„ I [cosh(2ax) - cos(2A;/x)l H{x), (29) 
1'^ + kj J 

\u{x)\^ = e^""" - ( ,A [asinh(2ox) ki^\u{2kix)\ H{x). (30) 
\a'^ + kf J 

Moreover, for the particular choice /cr = a we get a more compact expression for the new 
potential V2{x) than for a generic kji that would appear if we would substitute the u{x) and w{x) 
of equations ([26]) , (f27|) in equation (p8|) : 



V2ix) = -a5{x) + F(x)(^4a2e2'*^[2(a2 - k]) cos{2kix){a^ - kje^""") (31) 
-4akj sm{2kix){a^ + kje^'''') + Sa^kje^^'^ij/la^ - kje^'''' - 2a^e^''^ cos(2A;/x)]^ 



Let us remark that, for the general case characterized by equations ()26p and ()27p as well as 
the particular ones described by equations (f29]l - (f3T]) . the scattering states of are mapped 
into the corresponding ones of //2) and the same happens for the bound state. Thus, it turns out 
that the spectrum of li2 will be equal to Sp(//o) = {-S'o} U [0, oo), i.e., the complex second-order 
SUSY transformations which produce a real final potential are strictly isospectral. 



4.3 Real case 

Let us take now two seed solutions ui, U2 in the form given in equation Q, associated to the 
pair of real factorization energies e2 < (-x ^2 > Their explicit forms, and the corresponding 
derivatives, are given by: 

Ui{x) = 6^=^^ + Die-^-"" + 2ai sinh(A:ix)[A^(-x) - F(x)], (32) 
u^{x) = -km{x) + 2kie''''' [1 - diH{x) + DidiH{-x)] , « = 1, 2, 

where di = a/ki. Similarly as in the complex case, the calculation of the Wronskian w{x) = 
W{ui,U2) of the two involved Schrodinger seed solutions is once again cumbersome, but a con- 
venient compact expression reads: 

W{x) = {ki - k2)uiU2 + 2k2Uie^^'' [1 - d2H{x) + D2d2H{-x)] 
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Figure 3. Potential difference AV{x) as function of a; (blue lines), induced by a real second order SUSY 
transformation for a = 2, ^2 = 1, = 1/2, Di = —1/2, D2 = 1. Note that two new levels were created 
above Eq = —2, at the positions €2 = —1/2 and ei = —1/8 (gray horizontal lines). 



- 2kiU2e''^'' [1 - diH{x) + DiaiH{-x)] . 
By employing this equation, it is straightforward to calculate the new potential through: 



V2{x) = -a5{x) + 



w 



w 



W 
W 



-a6{x) + 



{kf - kl)uiU2 



w 



+ 



(kj - k^){uiU2 + U[U2) 



w 



Concerning the spectrum of H2, once again the scattering states of Hq are mapped into the 
corresponding ones of H2- As for the discrete part of the spectrum, several possibilities are 
worth of study. 

(i) Creating two new levels. Let us suppose first that ei ^ Eq ^ €2- In order that ti)(a;) 
do not have nodes, the two factorization energies must be placed either both below (for 
k2 > ki > a) or both above Eq (for ki < k2 < a). Moreover, according to the chosen 
ordering €2 < ei, the solution U2{x) must have one extra node with respect to ui{x) [37] . 
In the domain k2 > ki > a (e2 < ei < £"0) this can be achieved by taking D2 < and 
Di > while for ki < k2 < a {Eq < £2 < ei) it must be taken 1^2 > and Di < 0. With 
this choice of parameters, it turns out that the two eigenfunctions of H2 associated to ei 
and €2, V'ei^ oiU2/w and i})^^^ cc ui/w, are square-integrable. Thus, 

Sp(lf2) = {e2,ei}USp(//o), 

i.e., two new levels have been created for H2, either both below the ground state of Hq (for 
k2 > ki > a) or both above Eq (for ki < k2 < a). An illustration of the last situation is 
shown in Fig. [3l where we have plotted the potential difference AV{x) for a = 2, A;2 = Ij 
ki = 1/2, Di = —1/2, D2 = 1. As a result of the transformation, two new levels were 
created above the ground state energy of Hq at the positions £2 = —1/2 and ei = —1/8 
(see the gray horizontal lines at Fig. [3]) . 

(ii) Creating one new level. This case arises from the previous one for D2 — ?• 0. Now it turns 

(2) 

out that -0^2 is not square-integrable anymore, meaning that 



Sp(if2) = {ei}uSp(i^o) 



Thus, in order to generate H2 a new level has been created at ei, above Eq for ki < a and 
below it for ki > a. 
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(iii) Isospectral transformations. These can be achieved from case (i) for Di = D2 — ?■ 0, where 

(2) (2) 

both ipei and ipe2 cease to be square-integrable so that ej Sp{H2), i = 1,2. Hence, 
SpiH2) = Sp{Ho). 

(iv) Moving the level Eq. This procedure is obtained from case (i), e.g., by taking €2 = Eq, 
D2 — )• 0, U2{x) oc ipoix), and ui{x) as given in equation (j32]) with Di < 0, ei > Eq. 

(2) (2) 

With this choice it can be shown that oc ui/w is not square-integrable but ipei does, 
meaning that 

Sp(i72) = {ei}U[0,oo). 

In a way, the level Eq has been moved up to ei for generating H2. 

(v) Deleting the level Eq. This can be achieved as a limit of the previous case for Di — t- 0. 
Now it turns out that lim U2/W = 00, i.e., ei Sp(H2), and hence 

x—^O 

Sp{H2) = [0,00). 
5 Conclusions 

We have employed the first and second-order supersymmetric quantum mechanics for generating 
new potentials with modified spectra departing from the delta well potential. The first-order 
transformation allowed us to change just the ground state energy level, while the second-order 
transformations enlarged the possibilities of spectral control, including the option of manipula- 
ting the excited state levels. On the other hand, it is important to remember that the first-order 
transformations induced in the new potential a delta term with an opposite sign compared with 
the initial one (physically the delta term changed from attractive to repulsive). Meanwhile, the 
second-order transformations generated a delta term with exactly the same sign as the initial 
one (the attractive nature was preserved under the transformation). These physical differences 
should be taken into account in the determination of the most appropriate transformation for 
building a potential model. We can conclude that supersymmetric quantum mechanics is a po- 
werful mathematical tool, which is quite useful for implementing the spectral design in physics. 
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